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Abstract 

We study the effect of parton angular momentum on the twist-four correction to the left-right 
asymmetry in the electron-deuteron parity- violating deep inelastic scattering (PVDIS). We show- 
that this higher-twist correction is transparent to the dynamics of parton angular momentum 
needed to account for the Sivers and Boer-Mulders functions and spin-independent parton distri- 
bution functions. A sufficiently precise measurement of the PVDIS asymmetry may, thus, provide 
additional information about the parton dynamics responsible for nucleon spin. 



PACS numbers: 24.85.+p,13.60.Hb,11.80.-m,11.10.St 



I. INTRODUCTION 



As a complement to the studies at high-energy frontier, measurements at the intensity 
frontier (or precision frontier) provide powerful tools in the search for physics Beyond the 
Standard Model (BSM). Observables such as the muon anomalous magnetic moment are 
measured to very high precision, and experimental results are then compared with theoretical 
predictions. To the extent that the latter are sufficiently reliable, any possible deviation 
would point to BSM physics. Alternately, these experiments can provide new insights into 
the dynamics of the Standard Model. 

Electron-deuteron parity violating deep inelastic scattering (eD PVDIS) is an excellent 
example of this class of studies. Historically, it provided the first experimental measurement 
of weak mixing angle 9w [J]- Nowadays, with the prospect of the Jefferson Laboratory 
12-GeV upgrade and the use of a new spectrometer called SoLID, the left-right asymmetry 
of PVDIS can be measured with 0.5% precision over the kinematic range 0.3 < xb < 0.7 
[2]. With this level of precision, one will be able to probe or constrain an interesting set 
of BSM scenarios, such as a leptophobic Z' boson[3j and supersymmetryjl], as well as to 
study hadronic physics effects which are yet to be fully understood, such as charge symmetry 
violation (CVC) and higher-twist (HT). 

The effect of HT [5] is a potentially important Standard Model correction that originates 
from the interaction between partons. This correction in general scales as (Q 2 ) _( - r_2 ^ 2 , with 
the twist t > 2, so its effect is enhanced at low Q 2 . In the framework of the operator 
product expansion (OPE), the higher- twist correction can be expressed as a convolution 
of a high-energy and low-energy piece; the former (embodied in the Wilson coefficients) 
can be calculated using perturbative methods, whereas the latter involves hadronic matrix 
elements that require understanding of non-perturbative QCD. Studying the higher-twist 
correction may help us in probing correlations between the confined quarks and gluons 
inside the nucleon, so it is interesting to explore HT matrix elements within various model 
approaches. One advantage of eD PVDIS process is that the HT contribution to the leading 
term in the PV asymmetry (defined below) arises from a single operator matrix element and 
can, in principle, be separated kinematically from the subleading terms that have a more 
complicated HT structure. With this motivation in mind, several previous works [6HH] have 
been carried out to study the twist-four (i.e. r = 4) correction to the left-right asymmetry 
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of eD PVDIS. In what follows, we report on a study that follows-up these earlier works. 

The study of HT may also shed light on another important issue, namely, the spin 
structure of the nucleon. Nearly twenty-five years ago, the EMC collaboration [T5] performed 
a DIS experiment with longitudinally-polarized muons on a target of longitudinally-polarized 
protons, obtaining a value for the structure function Qi{xb) over the range 0.01 < xb < 0.7. 
After extrapolating to the low- and high-x^ region, the collaboration obtained a value for 
the leading moment of Qi(xb) that contradicted the Ellis- Jaffe sum rule [16] and implied 
that that the spin of proton is not built up entirely from the quark spin. The result has been 
confirmed by a variety of subsequent studies. A key question in nuclear physics research 
has, thus, become explaining in detail the source of nucleon spin in terms of QCD degrees 
of freedom. 

From a theoretical perspective, arriving at a decomposition of the nucleon spin in terms 
of gauge-invariant matrix elements of local operators that afford a straightforward partonic 
interpretation has been a vexing problem, and different approaches have been pursued over 
the years [TOUTi] . In each case, reference is usually made to the interpretation in the light- 
cone - gauge dependence notwithstanding - given its historical importance for thinking 
about parton dynamics. However, while the meaning of the quark helicity is gauge invariant, 
the relative importance of other aspects of partonic angular momentum (gluon helicity and 
quark and gluon orbital angular momentum) in general vary with the choice of gauge and 
even definition. Nonetheless, it is interesting to ask how different observables may probe 
different aspects of partonic angular momentum and to do so in a way that is both gauge- 
invariant and as insensitive as possible to a particular angular momentum decomposition. 

In this respect, we will study HT in the context of light-cone quantization. In early 
work within this framework, it has been shown that one particular component of parton 
angular momentum - identified as quark orbital angular momentum (OAM) under light-cone 
quantization using light-cone gauge - is responsible for the non-zero value of Sivers function 
and Boer-Mulders function [TTl HE] in semi-inclusive deep inelastic scattering (SIDIS) [19]. In 
light of these results, it is also interesting to study how the inclusion of the same component 
of parton angular momentum modifies the current model predictions for HT corrections to 
eD PVDIS. Indeed, in all the previous studies of eD PVDIS, only the Fock component of 
the nucleon wavefunction with zero parton OAM has been included. 

After including quark OAM in the light-cone amplitudes, we observe a rather non-intuitive 
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phenomenon: although the absolute magnitude of individual non-zero quark OAM contri- 
butions can be significant, they largely cancel against each other. We will argue that this 
cancelation is largely independent of the detailed model for the relevant light-cone ampli- 
tudes. As a result, the twist-four correction to PVDIS is almost transparent to the inclu- 
sion of quark OAM. In contrast, other hadronic quantities, such as the parton distribution 
functions (PDF), Sivers function, and Boer- Mulders function, manifest non-negligible de- 
pendence on quark OAM. Generalizing from the particular choice of light-cone quantization 
and light-cone gauge, we thus conclude that whatever features of parton angular momentum 
are responsible for the observed behavior of the PDFs, Sivers, and Boer-Mulders functions, 
they should have a relatively minor impact on the HT correction to eD PVDIS of interest 
here. Moreover, any deviation from the light-cone predictions obtained here and in previous 
works [BHS] - should they be observed expermentally - would signal the importance of other 
aspects of parton angular momentum and/or higher Fock space components of the nucleon 
wavefunction. 

The discussion of the computation leading to these observations is arranged in the fol- 
lowing order: in Section II we summarize the relevant results of the general formulation of 
the twist-four correction to eD PVDIS; in Section III we introduce the light-cone wavefunc- 
tion with quark OAM-dependence; in Section IV we present the analytic expressions of the 
hadronic matrix elements needed for the twist-4 correction, and demonstrate the generic 
cancelation between non-zero quark OAM components; in Section V we present the nu- 
merical results using one specific choice of nucleon wavefunction, and discuss their physical 
significance. Detailed formulae appear in the Appendix. 



II. HIGHER-TWIST IN PVDIS: GENERAL FORMULATION 

Here, we review the well-known results for the twist-four correction in eD PVDIS. We 
will simply quote the central equations that are relevant to our study without any derivation 
and refer the reader to Refs. [H [8] for the details. 

In eD PVDIS, longitudinally-polarized electron beams are incident on unpolarized 
deuteron targets. One measures the PV right-left asymmetry 

da R -da L 

Arl - —, — \y) 

aa R + cLctl 
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Figure 1: Kinematics of e-D PVDIS: a deuteron of momentum P interacts with an incoming 
electron of momentum k via an exchange of a single photon or Z-boson, and breaks into hadrons 
which are denoted collectively as X 

where dan/L is the differential cross-section for the scattering of the right/left-handed elec- 
trons. At the one-boson exchange (OBE) level, the leading parity-violating piece comes 
from the interference between photon and Z-boson exchange diagrams (see Fig |g . The low- 
energy Z-boson exchange interaction can be described by the following effective 4-fermion 
interaction: 

G F _ 

Cpv = ^|[e7 /i 75e(Ci u M7 At M + Ci d d^^d) + e^e{C 2u w-f^u + C^dj^d)) (2) 
where, at tree level, we have: 

C lu = -^ + ^sin 2 ^ (3) 

Cu=\- \&?0 W (4) 

C 2u = ~ + 2sin 2 w (5) 

C 2d = \ - 2sm 2 9 w (6) 

Neglecting contributions from sea quarks, assuming charge symmetry (v,y = d v , etc. with 
qy being the valence quark PDF of nucleon N), the leading-twist SM prediction is given by 
the Cahn-Gilman formula[9]: 

Arl = \ [(2C lu - C u ) + (2C 2u - C 2d ) (7) 

2V27ra5 1 + (1 - yy 

where Q 2 = —q 2 and y = P ■ q/P ■ k. 



To include corrections from possible BSM and as well as other SM pieces, we can 
reparametrize the Cahn-Gilman formula [7]: 

G F Q 2 3,. l-(l-y) 2 , . 

ARL = -W^aS [ai + " 2 i + V-y)> ] (8) 

with dj = — (2Ci U — Cid)(l + Ri). Here, Ri describes any deviation of the Cj from the 
expressions in Eqs. ([3]) to (j6]), including both SM and BSM corrections. In this paper we 
concentrate on Ri T , namely the higher- twist correction to a\. 

Bjorken and Wolfenstein (20J, EJJ showed that, if one assumes isospin symmetry and 
neglects sea quark contributions, then there is only one matrix element that contributes to 
Ri T (for a detailed review of these arguments in a more modern context, see Ref. [7J). This 
observation significantly simplifies the theoretical interpretation of the asymmetry, allowing 
us to concentrate on one particular matrix element without needing to to disentangle the 
contributions from many different sources. In brief, the Bjorken and Wolfenstein argument 
works as follows: Arl arises from the interference between the electromagnetic and weak 
neutral currents. First, one can decompose both currents into an isoscalar S and an isovector 
V term. The matrix elements of the S x V cross-term vanishes because deuteron is an 
isosinglet. Furthermore, at leading twist, we have (SS) = (VV). Therefore, the difference 
between (SS) and (VV) that enters hadronic tensor W^ u 

Kd(P,l) = 7t4t- I d ' zeiq ' Z (D(P)\u(z)Yu(z)d(0)rd(0) + (u^d) \D(P)) (9) 
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with Mjj being the mass of deuteron, is the only matrix element giving a HT correction to 

Below, we will compute the matrix element <Efy using an expansion of string operators 
[22] in order to extract the twist-four piece; the latter is expressed in terms of the deuteron 
twist-four distribution function Qd(xb), which will be computed in Section IV. 



III. THE LIGHT-CONE AMPLITUDES 

The main challenge in proceeding from ^ is our ignorance of the details of the nucleon 
wavefunctions. As QCD is non-perturbative at the hadronic scale, analytical expressions for 
the wavefunctions are unknown. At present, lattice QCD can provide only HT contributions 
to structure function moments and not the x^-dependence of the R± T that is of interest 
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to the SoLID experiment. Consequently, one must turn to various models that seek to 
incorporate non-perturbative dynamics. Previous works on R^ T include the use of MIT bag 
model [7] and isotropic light-cone wavefunctions that contain both quark and gluon Fock 
components [8]; their results yield similar shapes for the rcg-dependence but differ somewhat 
in magnitude, with a maximum i?f T of 0.003 ~ 0.005 at 0.2 < x B < 0.7 for Q 2 = 4GeV 2 , 
which is a little bit lower than the achievable precision level in the SoLID experiment. 

In this work we study how the inclusion of additional parton angular momentum might 
modify the R± T prediction. For this purpose, we adopt the formalism developed in Ref. |23j, 
starting from a light-cone formulation of quark states which is equivalent to the well-known 
"infinite momentum frame" point of view that gives the PDF its intuitive meaning as a 
parton momentum probability distribution [5]. We then perform a light-cone expansion of 
the nucleon state, retaining only the portion of Fock space containing three valence quarks 
with all possible quark OAM. To illustrate, we consider a spin-up proton. Its three valence 
quarks can form a total helicity of ±1/2, ±3/2; therefore in order to keep the total proton 
spin in z-direction to be 1/2 we need to assign different z-component quark OAM (i.e. l z ) 
for different combinations. 

A spin-up proton state, then, can be parametrized as the follows: 



\p f) = \p t) iz -° + \p I)'*- 1 + \p t) /z__1 + \p t) lz_z (10) 

,abc r 

\Pt) h= ° = -^y[M s ](V' (1) (l,2,3)+i(^-*i'A5)^ 2) (l,2,3))x 

«!t(i){«k(2)4i>(3) - <( 2 X( 3 )} 1°) ("J 

,abc r 

|P|)' Z=1 = ^y[M 3 ](fc 1 W 3) (l ) 2,3) + ^ {4) (l ) 2,3))x 

« t (l)<(2)4;(3) - 4t(l)<(2)<(3)) |0) (12) 

,abc r 

\PA) h= ~ l = ^y[ J DA 3 ](-^ 5 )(l,2,3))(ni t (l>4 t (2)< 1 .(3)x 

-ti+ t (l)4 t (2)ui f (3) |0» (13) 

l z =2 e 

«L (2)4(3)) |0) (14) 



with 



r abc 



\P^t= 2 = -j= /[M 3 ]Ai^ W (l,2 > 3)(tt^(l)4(2)ttt t (3)x 



where kf ± = kf ± ik^, while u ai (l) means the creation operator of an up-quark (same 
for down-quark) with color a, spin % and momentum k\ etc, satisfying the light-cone anti- 
commutation relation: 

{uai(p), = 2p + (2n) 3 5 ab S lJ 5(p + - p /+ )6^(p ± - ± ) (15) 

The integration measure is 



r- dx 1 dx 2 dx 3 d 2 k 1± d 2 k 2± d 2 k 3± 

[DX:>1 = ^ vafes; — Wf — 27tS<1 - 11 - 12 - x 

(2,r) 2 * 2 (fc 1± + k 2 ± + h±) (16) 
The proton wavefunction amplitudes {^^...tp^} are generally unknown functions. Al- 



though the expansion (11)~(14) is generic, the explicit form of if)® is model-dependent. In 
this work, we chose the form of if)® derived in Ref. [T7| by starting from the static solution 
of a constituent quark model [21] (which works well in predicting many electroweak prop- 
erties of the baryons) and applying a Melosh rotation to the solution to obtain non-zero l z 
components [25]. This choice of proton wavefunction is used to predict the first moment of 
Sivers function, and turns out to agree fairly well with the experimental measurements from 
HERMES and COMPASS 



IV. MATRIX ELEMENTS BETWEEN NUCLEON STATES 

Following [22], in order to obtain the twist-four distribution function Qd{x) we need to 
evaluate the matrix elements between state \D(P)) of the following operators: 

Q A (b,z) = : u(b lZ )t a ^ 5 u(b2z)d(b 3 z)t a ^d(b 4 z) : 

Q v (b,z) = ■.u{b 1 z)t a ^u{b 2 z)d{b 3 z)t a ^d{b 4 z) : (17) 

where z is a coordinate on light cone, and the parameters b = {&i, b 2 , b 3 , 6 4 } characterize the 
light-cone separation between quark field operators. 



When computing the matrix elements of Qva i n Eq. ( 17) we assume an incoherent impulse 



approximation in which the incoming photon strikes only one of the two nucleons; hence, 



matrix elements of the operators (17) can be related to the same matrix elements taken 



between proton states (or equivalently between neutron states, given isospin symmetry). 
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Also, since the quantities we compute do not depend on the proton spin, we can take it to 
be +1/2 along the z-direction without loss of generality. 



Now, starting from the operators (17), we define two distribution functions Q±(x^) via 



{Pip) tl {Qv(b, z) ± Q A (b,z)} \P(p) t> = (p ■ zf ! f[ dx^5(J2^)e-^ z ^ b ^Q ± (x^ 

J k=l i 

(18) 

with x^ collectively representing {x^, xg 2 , x^ 3 ,x^}, the light-cone momentum fractions: £ 4 + = 
x^p + . Meanwhile \P(p) T) is the spin- up proton state with momentum p. Substituting 
( 1 1 )~( 14 ) into (18) we are able to express Q±{x^) in terms of the proton wavefunction 
amplitudes. It is easy to observe that only diagonal terms, (i.e. terms with the same l z 
in initial and final states), could give non- vanishing contributions. After a rather lengthy 



derivation with the aid of Eq. (A2), we obtain 



327T 3 f d 2 £i± d 2 C,4± 
Q±(xt) = g- J j^---j^S(-x^)6(x^)6(-x^)e(x u )9(l-x^-x u )x 

5 2 (Zi± + £ 2 ± + £ll + Ol) E ^(-6, -6, 6, &) (19) 

where the explicit formulas of ^ are given in Appendix |b| 

The proton twist-four distribution function can now be expressed in terms of the Q± 
(refer to Eq. (42) of Ref. [S] after some rearrangement): 

/i tt 4 ^ x 
^ — — -<5(y^x ? J{(x 6 +x 6 )<5(xfl + ar ft +x 6 ) - x 6 <5(x B + x a ) 

-afctffo, - x B )}[(l + PuP2s)Q+(xt) - (Pn + P 3 4)<5-(^)] (20) 

Here Pij is the permutation operator, e.g. Pi 2 Q+(x^, x&, x^ 3 , a^J = Q + (x^ 2 , x^, x^ 3 , x^). 
The deuteron twist-four distribution function Qd(%b) can be expressed in terms of Q p (xb) 
through an incoherent impulse approximation [27], which says that a general deuteron 
hadronic tensor can be related to the corresponding hadronic tensors of proton and neu- 
tron by: 

M D W^(p, q) « M N Wf& q) + M N Wr(^, q) (21) 



where is the mass of nucleon. In the equation above each hadronic tensor is multiplied by 
the particle's mass, because following Eq. ^ the hadronic tensor we defined has dimension 



-1. Now we can express both sides of Eq. (21 ) in terms of dimensionless structure functions 
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{Fi(x B )}. Using isospin symmetry and the fact that Q(x B ) is proportional to x B Fi d (x B ) 
(see Eq. (34) of Ref. [HI), we obtain [53]: 

(22) 



^Qd{xb/2) ~ + Q n {x B ) « 2Q p (a; B ) 



Finally, following the logic of Ref. [8], one can the derive the twist-four contribution to 



Rr. 



xQ D {x B ) 



Q 2 5(1 - f sm 2 6 w ) u D (x B ) + d D (x B ) 



(23) 



with q B {x B ) being the parton distribution function for quark of flavor q in the deuteron 

(D(P)\ q(z)tq(-z) \D(P)) = 2(P ■ z) j dxe 2l{p - z)x q D {x) (24) 

Note that we neglect the logarithmic Q 2 -dependence of the structure functions in this anal- 
ysis. We can express q B in terms of PDF of the proton and neutron again by the impulse 



approximation (21), but now comparing the structure function F 2 (x B ) on both sides, which 



is proportional to x B 1 q(x B ). The result is: 



q D (x B /2) « q p (x B ) + q n (x B ) 



(25) 



where q p (x) and q n (x) are defined as in Eq. (24) but for proton/neutron states. Furthermore, 
neglecting CSV effects we have: 



Uu(xb) = d p (x B ),d n (x B ) = u p (x B ) 



(26) 



Therefore, it is sufficient to just calculate u p (x B ) and d p (x B ) using the proton light-cone 



wavefunction (11)~(14). Using (A3) and (A4) , we can compute the quark PDFs of the 



(spin- up) nucleons by calculating the matrix element on LHS of Eq. (24) with nucleon 



states, and compare it with the form on RHS to extract the PDFs. Same with the twist-four 
distribution functions, only terms diagonal to l z survive, so we can separate the result into 
components of different l z as the following: 

u p {x B ) + d p (x B ) = d n (x B ) + u n (x B ) 

= J^yJ o dx ± j d 2 k 1± d 2 q 1 _e{l-x B -x 1 )Y,^{q,l,2) (27) 



where the functions A z (q, 1, 2) are given in Appendix \B 
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/ 1 1 1 ] finnfn Knto y"y~i nQr nr T/~» 

fc 2 — UjH^X LUIlLIlUULfcJ IIlUoLiy LU 


anrl lpcc ^~ <■ \ /O - ^ - ■ urVi i 1 o / — 1 /"» n 1 t r r'i~»"n'4"T"iT^n'l"oc ^ /"i "T" 
dllU. Icoo LU , Willie 6^ — J- ^■^■'"'■J' ^UllLllUUlCo LU . 




^ Contribution to Qp{xs) 


Contribution to Qp(x B ) 


Subdominant 


Dominant 


+1 Subdominant 


Dominant 


-1 All 


Zero 



We now proceed to show that a partial cancelation occurs between contributions of l z 



+1 and l z = — 1. For this purpose, we combine (19) and (20), together with the fact that 
4>i (q, I, q' 1 1')* = ipi (q 1 ■, I', q, I), to simplify the expression of Q p (x B ) as: 



Q P (x B ) = Qp{x B ) + Q p (x B ) 



(28) 



where 

Q^b) 



3 pi 4 



I IT dx ^ x ^ - x 6 + ^3 - x ^) e ( l - x & - a; 54 ){ 

5(x B -x fl ) <5(xB-a; 54 ) : 5(x B + x g3 

H } x 



5(x B - Zft + 

x^ 2 x^ 3 



X & ( X & X £,3 ) X £,3 ( X & ^£3 ) 



x £i x U 



I n C% 52 (^i - £u + &3 - a 4 ) £ Re^+(6, 6, 6, &) 

^ i=l ^ ^ «z 



(29) 



Q p ( x b) = / JJc%5(:r 5l -x i2 -x& + x u )0{\-x i2 -x^){ 5 ^' 1 ' 



8(x B — x^ 4 ) S(x B — XfJ <5(xs — + 



-}x 

11 T^^ 1 " ^ 2 " + f±4) £ R<3 ^ 6, ( 30 ) 



i=l 



First we qualitatively analyze the contribution from each /^-component to Q^(x B ). This 
can be done by simply referring to Eqs. (B1)~(B8) of the Appendix B. The result is sum- 
marized in Table [ij We observe that the l z = +1 (-1) piece contributes mainly to Q~ (Qp)- 
Also notice that we do not include the l z — 2 component as its effect is tiny. 

Next we study the rough behavior of Qp (x B ) with respect to x B . For this purpose let us 
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Figure 2: The rough feature of Qp(x) ignoring the details of nucleon wavefunction amplitudes. 



The constant C is defined in Eq. (31). 
naively take: 



/ TT 7^?< j2 (-) Re $ t « constant = C (31) 
J 7" i ( 27r ) 



ft 

1 (27T) 

assuming the function above is well-behaved with respect to This approximation 

simply means that we do not care about the details of the proton wavefunction amplitudes. 



Under this approximation, the numerical integration (29) and (30) can be performed quite 
trivially, and the result is shown in Fig 2 One can see that a sign change of Qp(xs) 
from negative to positive occurs near xb = 0.4, implying that (xb) and Q~(xb) will 
have different signs for xb > 0.4. Therefore, according to Table [I], the contribution to 
Q p (xB)irom l z — 1 and l z = —1 should partially cancel other in this region. Furthermore, 
since the argument above does not depend on the details of nucleon wavefunction (as long 
as it is well-behaved), this feature of partial cancelation should be generic. The numerical 
results we present in the next section that include the full-dependence of t/>W on the £j, 
indeed, exhibit this near cancellation. 

V. NUMERICAL RESULTS AND DISCUSSION 



Eqs. (29) and (30) are our starting point for the numerical evaluation of Q p (xb), which 
involves an eight-fold integration. To perform this integration, we adopt the Monte Carlo 
numerical integration called Divonne contained in the CUBA Library, which is an algorithm 
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Figure 3: (color online)The Twist-4 correction to R\ at Q 2 = 4GeV 2 . The blue dashed curve shows 
the l z = contribution; purple dot-dashed curve shows the l z = 1 contribution; brown dot-dashed 
curve shows the l z = —1 contribution; the red solid curve is the sum of all. l z = 2 contribution is 
negligible and therefore not included. 

package designed for multi-dimensional numerical integration [2S]. For each l z component, 
we evaluate the value of Q p (xb) at a series of discrete {xb,%}, and then link them together 
using a best-fit line. Also, we take a s = 0.5 at lGeV following the renormalization group 
(RG) prediction of the running coupling constant at 4-loop order together with a 3-loop 
threshold matching, with the quark thresholds taken to be M c = 1.5 GeV and M& = 4.7 
GeV respectively [2"§] . 

Our main result is shown in Fig. [3J which gives versus x' B = 2x B at Q 2 = 4GeV 2 . 
First, let us compare this outcome with that of Refs. [7] and [SJ. It turns out that all three 
calculations predict similar curve shape for Ri T , only with slightly different positions of 
peak and zero-point. Concerning the magnitude, our work predicts a maximum absolute 
value |-R^ T | ~ 2.6 x 1CT 3 between 0.2 < x' B < 0.7, which is smallest in magnitude among 
all the three predictions, and is about a half of the size to that of Ref. [5j. This is un- 
derstandable because the authors include a 3-quark+l-gluon Fock-space component whose 
contribution is comparable in magnitude to that of the pure 3-quark state. Nonetheless, all 
three calculations suggest that |-R^ T | lies below that of the expected SoLID precision. 

Next we study the OAM-dependence in detail. To that end, we first introduce some 
nomenclature: in the following, we will use the notation (|/ 2 | ® \l z '\), which denotes a generic 
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Figure 4: (color online)The unnormalized QDF of spin-up proton, splitted into contributions from 
different l z components. Blue thick-dashed curve shows contribution from l z = component; 
purple dot-dashed curve shows contribution from l z = 1 component; brown dot-dashed curve 
shows contribution from l z = — 1 component; green thin-dashed curve shows contribution from 
l z = 2 component; red solid curve is the sum of all contributions. 

matrix element taking between two hadronic states, of which one of them has absolute value 
of quark OAM in z-direction being \l z \ and the other being \l z i\. 



From our arguments at the end of Section IV, we expect that although l z = ±1 individu- 



ally contribute a significant amount to Q p (xb), they should largely cancel against each other 
for xb > 0.4, making the total (1 £g> 1) contribution rather small, and therefore leaving the 
(0 Cg> 0) contribution as the dominant piece. This expectation is born out by the curves in 
Fig. [3] The purple dot-dashed curve and brown dot-dashed curve curves give the individual 
{l z — 1) <g) (l z = 1) and {l z = —1) Cg) [l z = —1) contributions, respectively , which exhibit 
the expected cancellation for x' B > 0.4. The blue dashed curve and red solid curve give the 
(0 Cg> 0) and total contributions, respectively. It is clear that the former dominates the total. 
This (0 £g> 0) dominance is a rather unique feature of the particular twist-four contribution 
of interest here, and one that is not shared by other diagonal matrix elements. For example, 
if one calculate proton quark PDFs (leading twist) using the same set of wavefunctions, the 
(0 ® 0) and (1 (g> 1) contributions are comparable; moreover, since they have the same sign, 
the two \l z \ = 1 pieces do not cancel each other (see Figj4]). 

On the other hand, we also note that there are hadronic matrix elements that depend 
crucially on the existence of non-zero quark OAM in light cone quantization. In partic- 
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ular, in Ref. [T7], the authors studied the Sivers function [30] and Boer-Mulders function 
[3"T] . which are examples of transverse momentum dependent parton distribution functions 
(TMDs), appearing in semi-inclusive deep inelastic scattering. Importantly, both distribu- 
tion functions depend on off-diagonal matrix elements of l z : the Sivers function is sensitive 
to (0 ® 1) while Boer-Mulders function is sensitive to both (0 ® 1) and (1 £g> 2). Simply 
speaking, the existence of non-zeo quark OAM is responsible for the non-vanishing values 
of the Sivers and Boer-Mulders functions. Combining this observation with our analysis of 
the HT matrix element, we conclude that the twist-four correction to eD PVDIS is essen- 
tially transparent to the parton angular momentum dynamics that generate the Sivers and 
Boer-Mulders functions. 

It is also interesting to study the impact of sea-parton dynamics on the behavior of the 
HT matrix element. To that end, we performed a qualitative analysis of the contribution 
made by the Fock space component containing 3 quarks + 1 gluon, using the general form 
suggested in Ref. [32] that includes non-zero gluon OAM. The authors of Ref. [S] computed 
the contribution of the 3q+lg state with l z = , which turns out to have a very similar 
shape with the l g = contribution of 3q-state. To our knowledge, however, there exist no 
explicit functional forms for the 3q+lg nucleon wavefunction with non-zero parton OAM. 
Consequently, our analysis is purely analytic at this point. Nonetheless, in contrast to the 
quark OAM case, we find no obvious pattern of cancelations between components of non-zero 
l z ■ 

Combining observations, we may draw the following conclusion: if a future eD PVDIS 
measurement yields a sufficiently precise determination of R^ T as a function of x' B , one 
can compare the experimental curve with our current theoretical prediction. A significant 
deviation from the predicted curve (e.g., the peak and zero-point are shifted by a considerable 
amount), could signal the importance of parton angular momentum dynamics beyond those 
responsible for the Sivers, Boer-Mulders, and spin-independent parton distribution functions. 

VI. SUMMARY 

The next generation of parity-violating electron scattering experiments are poised to 
probe both possible BSM physics as well as novel features of hadron and nuclear struc- 
ture. In this work, we have studied one particular hadronic effect, namely, the twist-four 
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Table II: The dependence on different quark light-cone OAM components of various distribution 
functions. 



Distribution Functions Dominant Contributor(s) Subdominant Contributor(s) 
Quark Distribution Functions (0x0), (1(8)1) (2(8)2) 
PVDIS Twist-Four Correction (000) (1®1), (2(82) 

Sivers Function (0(81) (1(82) 

Boer-Mulders Function (0(81), (1(82) 



contribution to 5i, the y-independent term in the PV asymmetry. Using a set of proton 
light-cone wavefunctions with non-zero quark orbital angular momentum, we evaluated the 
twist-four contribution as a function of xb, identifying the contributions from different OAM- 
components. Our total for the correction R^ T is similar in both shape and magnitude to 
those obtained in previous works, indicating that higher-precision than expected with the 
SoLID experiment would be needed to discern this HT effect. An effort to achieve such 
precision may be worthwhile, because R^ T appears to be rather unique, in the sense that 
it is not significantly affected by the parton angular momentum physics responsible for the 
existence of some other DIS observables such as the Sivers and Boer-Mulders functions. 
Thus, by combining the results of a more precise measurement of the asymmetry with mea- 
surements of other distribution functions, it is possible to probe complementary aspects of 
parton angular momentum and, perhaps, shed new light on the role of angular momentum 
in the structure of the nucleon. 
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Appendix A: Matrix Elements of Two and Four-Fermion Operators 

In this section we present matrix elements of two-fermion operators {v)u and d)d) and 
four-fermion operators {v)ud)d) between nucleon states. For this purpose let us consider 
two arbitrary components of proton light-cone wavefunction defined as the following: 

f abc r 

= ^/[ M 3]V'a(l,2,3)^ Ai (l)^ A2 (2)4 A3 (3)|0) 

f abc r 

|^> = -j=J [M 3 ]^(1,2,3)^ A ,(1)^ A ,(2)4,(3) |0) (Al) 

It is straightforward to work out the matrix elements of the four-fermion operator between 
these two states: 



{i>a\Ap{q) U i'fs{d)d] x {l)dj'X>{l') |#) = —^-(SwSjji - Sij'Si'^SxsxSx^X' sj XqXiX'gX't 

5(x q + xi- x' q - x'^S 2 ^ +l±- q' ± - l'±) J dx!d 2 k 1± 6(l - x 1 - x q - Xi)S 2 (k 1± + q± + l±) 
(Sx lP Sx 2 x' 2 Sp'x' 1 ^a(^ !> WpW, !> l ') + h 1 \' 2 Sx 2P S p 'x' 1 ^*a(^, q, l)i>p{q , 1, 1') 
+Kph 2 x> 1 Sp>x> 2 Tp* a (q, 1, 0^(1, <?', + Sxx^pS^rfcil, q, Z)^(l, g', /')) (A2) 

and those for two-fermion operators: 

(V'al djx( l ) d j'X'( 1 ') \4>p) = ^Xi5(xi-x' l )5 2 (l±-l 1 ± )5x :} x5x' :t x'5jj' J dx 1 dx 2 d 2 k 1 ±d 2 k 2 ± 

S(l-xi-x 2 - xi)S 2 (k 1± + k 2± + l±)(5x 1 x> 1 5x 2 \ 2 ipi3(^, 2, 1) 
+5 Xl x> 2 h 2 x>M 2 i 1,0)^(1,2,0 (A3) 

(V'al u] x {l)u r x'{l') \4>p) = ^xi5(xi - x^S 2 ^^ -f ± )6x 3 x> 3 8jj> J dx x dx 2 d 2 k xl _d 2 k 21 _ 

5(1 - X! - x 2 - ^)5 2 (^i± + k 2 ± + f±)(5 Xl xS\ 2 \ 2 Sx\> 1 '4>*a(l, 1, 2)^(i, 1, 2) 
+5a 1 a^a 2 a5va' 1 C(1 ! 2)^(Z, 1, 2) + (WW^'Ca 1, 2)^(1, /, 2) 
+5a 1 a'/a 2 a5a^C(1> 2)^(1, 1, 2)) (A4) 

Appendix B: Complete formulae for quark PDFs and Q p {xb) in terms of proton 
wavefunction amplitudes 

In this section we present explicit expressions needed to compute the quark PDFs and 
the twist-four distribution function. 
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The distribution functions Q^(xb) in Eqs. (29) and (30) are expressed in terms of 
4>^(q, I, q', I'), which have the following expressions: 

rl>t=o(Q,l>M = 2^ 2 >(q,l,l)^ 2 \q',l,l') (Bl) 
tf£=o(?, 4, n = 2{V (1 ' 2) *(1, q, 0^ (M) (1, + ^ (1 ' 2) *(9, «, 1)V> (1 ' 2) (9', J', 1) 

+^ (1 > 2) *(i, z, ?)v (1 ' 2) (? / , r, 1) + v (1,2) *(?, 1, i)^ (1 ' 2) (i> ?') 

+^ 1 ' 2) *(l,/,g)^ (1 ' 2) (l, (B2) 

^=i(g,W,0 = 2^ (3 ' 4) *(i,g,/)^ (3 ' 4) (i,g / ,n (B3) 
V£ =1 (g, J, </, = 2{^>(q, 1, Z)^ (3 ' 4 V, 1, + V> (3 ' 4) *(Z, ff, 1)^ (3 ' 4) (^ 9, 1) 
+^>(l, 1, ? )V (3 ' 4) (r, q\ 1) + ^ (3 ' 4) *(Z, q, 1)^(1', 1, q') 
+^ 3 ' 4 >(l,l,q)^\l',l,q')} (B4) 
^=-i(g, l, q', = 2{^( 5 ' 5 )*(g, 1, Z^ {5 ' 5) (</, 1, + ^ (5 ' 5) *(1, ff, 0^ (5 ' 5 V, 1, 

+^ (5 ' 5) *(g, 1, /)^ (5 ' 5) (1, g', + V> (5 ' 5) *(1, q, /)^ {5 ' 5) (1, 0} (B5) 
= (B6) 

z ' ^ = 2{^ (6 ' 6) *(g, 1, o^ (6 ' 6) (?', 1, + ^ <M> *(i, 9, ov (e ' B V. 1, n 

+^ e >(q, 1, ity (6,6) (l, q', I') + V> (6 ' 6) *(1, q, l)^ 6 \l, q', l')} (B7) 

^(g, /,(/,/') = (B8) 

The symbol "1" denotes the four momentum k% = {xip + ,ki±) which is given by x% = 
1 — x q — xi = 1 — x' q — x\ and fcu_ = — <f_i_ — l± = — <f|_ — The definitions of ijj^> are the 
following: 

^(1,2,3) = ^(1,2,3) +i(kfk v 2 -kfk^ 2 \l,2,3) 
^ (3 ' 4) (1,2,3) = k+ ± ^\l,2,3)+k+ ± ^ 4 \l,2,3) 
^ (5 ' 5) (1,2,3) = -k- ± ^ 5 \l,2,3) + k~ ± ^ 5 \l,3,2) 

^(1,2,3) = fc 1 \(^( 6 )(l,3,2)-A; 4 ± ^ 6 )(l,2 ) 3)) (B9) 



On the other hand, the quark distribution functions in (27) are given in terms of 
A lz (q, 1, 2), which look like the following: 
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A l *=°(q, 1, 2) = V (1,2) *(9, 1, 2)^ (1 ' 2) (g, 1, 2) + 2^(1, g, 2)^ (1 ' 2) (1, g, 2) 
+^ 2 >(q, 2, l)^ (1 ' 2) (g, 2, 1) + V (1 ' 2) *(l, 2, g)V (1,2) (?, 2, 1) 
+^ {1 ' 2) *(g, 2, 1)V (1 ' 2) (1, 2, q) + 2^ (1 ' 2) *(1, 2, g)^ (1 ' 2) (l, 2, g) 
+^ {1 ' 2) *(l,g,2)^ 1 ' 2) (2,g,l) (BIO) 

^ =1 (g, 1, 2) = 2^ 3 > 4 >*(g, 1, 2)^ 3 > 4 )(g, 1, 2) + ^ 3 ' 4 )*(1, q, 2)^(1, g, 2) 
+^ 3 ' 4 )*(2, g, 1)^ 3 < 4 )(2, g, 1) + ^ 3 < 4 )*(2, 1, g)^ 3 ' 4 )(2, g, 1) 
+^ 3 < 4 )*(2, g, 1)^(2, 1, g) + ^ 3 < 4 )*(2, 1, g)^ 3 ' 4 )(2, 1, g) 
+^ 3 ' 4 )*(1, 2, g)^ 3 ' 4 >(l, 2, g) + V (3 ' 4) *(g, 1, 2)^ 3 ' 4 )(g, 2, 1) (Bll) 

^=- 4 (g, 1, 2) = ^ (5 ' 5) *(g, 1, 2)^ 5 ' 5 >(g, 1, 2) + V (5 ' 5) *(l, g, 2)^(g, 1, 2) 
+^ (5 ' 5) *(g, 1, 2)^ (5,5) (1, 2) + ^ (5,5) *(1, q, 2)^ (5 ' 5) (1, 4, 2) 
+^ (5 ' 5) *(1, 2, g)^ (5 ' 5) (l, 2, g) + ^ (5,5) *(1, 2, g)^ (5 ' 5) (2, 1, g) (B12) 

^= 2 (g, 1, 2) = V (6 ' 6) *(g, 1, 2)^(g, 1, 2) + ^ (6 ' 6) *(1, q, 2)^\q, 1, 2) 
+^ (6 ' 6) *(g, 1, 2)^(1, g, 2) + ^(1, g, 2)^ (6 ' 6) (1, q, 2) 
+V> (6 ' 6) *(1, 2, g)^ (6 ' 6) (l, 2, g) + ^ 6 ' 6 >*(1, 2, g)^ (6 ' 6) (2, 1, g) (B13) 

with g = (xbP + , gi), ^2 = 1 — %b — x\ and &2± = — gj_ — ki±. 
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